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Abstract 


The  hydrodynamic  forces  on  the  Dorado  semi-submersible  remote  minehunting  vehicle  are 
estimated  and  provided  in  a  form  that  can  be  used  in  a  six  degree-of- freedom  vehicle 
simulation.  The  forces  apply  to  either  old  or  new  versions  of  the  control  planes.  The  lift 
(from  flow  curvature  in  a  turn)  and  drag  of  the  self-aligning  fairings  on  the  snorkel  mast  are 
accounted  for  as  a  function  of  vehicle  depth.  Hydrodynamic  coefficients  are  used  for  all  but 
the  mast  effects,  where  a  rational  flow  model  is  used.  Deficiencies  in  the  control  derivatives 
estimated  by  the  DERIVS/DSSP20  programs  are  identified  and  overcome. 


Resume 


Les  forces  hydrodynamiques  exercees  sur  le  vehicule  semi-submersible  Dorado  utilise  pour  le 
deminage  a  distance  sont  evaluees  et  fournies  sous  une  forme  utilisable  dans  une  simulation 
de  vehicule  a  six  degres  de  liberte.  Les  forces  sont  appliquees  soit  a  la  nouvelle  soit  a 
l’ancienne  version  des  plans  de  commande.  La  poussee  (due  a  la  courbure  de  l’ecoulement 
dans  un  virage)  et  la  trainee  du  carenage  a  auto-alignement  sur  le  mat  du  schnorkel  sont  prises 
en  compte  sous  la  forme  d’une  fonction  de  la  profondeur  du  vehicule.  Des  coefficients 
hydrodynamiques  sont  utilises  pour  tous  les  effets,  sauf  celui  du  mat,  pour  lequel  un  modele 
d’ecoulement  rationnel  est  utilise.  Des  lacunes  concernant  les  parametres  de  commande 
evalues  par  les  programmes  DERIVS/DSSP20  sont  identifiees  et  corrigees. 
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Executive  summary 


Introduction 

DRDC  Atlantic  has  developed,  and  is  now  supporting  a  Technology  Demonstration  Project  of 
a  Remote  Minehunting  System  (RMS).  The  RMS  uses  a  remote  drone  (called  Dorado)  to  tow 
a  sonar  at  depths  up  to  200  m  for  minehunting  and  route  surveying  purposes.  DRDC  Atlantic 
has  also  developed  a  six  degree-of-freedom  simulator  (SIMRMS)  for  Dorado,  its  tow  cable, 
and  its  towfish.  The  simulator  requires  the  hydrodynamic  characteristics  of  the  RMS 
components  as  input.  This  report  estimates  these  characteristics  for  the  Dorado  vehicle. 


Principle  Results 

The  Dorado  hydrodynamic  characteristics  are  estimated  as  coefficients  and  rational  flow 
models.  Propulsion  effects  are  not  dealt  with  in  this  report  since  further  changes  to  the 
propulsion  system  are  expected.  During  the  calculations,  deficiencies  were  identified  in 
programs  which  have  been  used  in  the  past  to  estimate  control  derivatives  (depth,  pitch,  yaw, 
and  roll  control  authority).  Improved  algorithms  which  correct  these  deficiencies  are 
presented. 

Significance  of  Results 

The  new  Dorado  hydrodynamic  model  should  allow  SIMRMS  to  function  with  improved 
accuracy.  SIMRMS  is  an  important  simulation  tool  in  planning  and  carrying  out  remote 
minehunting  and  route  surveying  operations. 

Future  Plans 

It  is  desirable  to  validate  the  improved  control  derivative  algorithms  through  full  scale  trials 
with  Dorado.  Such  trials  are  planned  but,  as  yet,  unfunded. 
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Sommaire 


Introduction 

R  &  D  pour  la  defense  Canada  -  Atlantique  a  mis  au  point  un  systeme  de  deminage  a  distance 
(SDD)  dans  le  cadre  d’un  projet  de  demonstration  de  technologie  et  il  fournit  actuellement  le 
soutien  a  ce  systeme.  Le  SDD  utilise  un  engin  telecommande  (appele  Dorado)  pour 
remorquer  un  sonar  a  des  profondeurs  pouvant  atteindre  200  m  a  des  fins  de  deminage  et  de 
leve  des  routes  de  navigation.  R  &  D  pour  la  defense  Canada  -  Atlantique  a  egalement  mis  au 
point  un  simulateur  a  six  degres  de  liberte  (SIMRMS)  pour  le  Dorado,  son  cable  de 
remorquage  et  son  sonar  remorque.  II  faut  fournir  au  simulateur  les  caracteristiques 
hydrodynamiques  des  elements  du  SDD.  Dans  le  present  rapport,  on  evalue  ces 
caracteristiques  pour  le  vehicule  Dorado. 


Principaux  resultats 

Les  caracteristiques  hydrodynamiques  du  Dorado  sont  evaluees  sous  forme  de  coefficients  et 
de  modeles  d’ecoulement  rationnel.  Le  rapport  ne  traite  pas  des  effets  de  la  propulsion,  etant 
donne  que  d’autres  changements  sont  prevus  au  systeme  de  propulsion.  Lors  des  calculs,  on  a 
identifie  des  lacunes  dans  les  programmes  utilises  par  le  passe  pour  evaluer  les  parametres  de 
commande  (profondeur,  tangage,  lacet  et  roulis).  Des  algorithmes  ameliores  qui  corrigent  ces 
lacunes  sont  presentes. 

Signification  des  resultats 

Le  nouveau  modele  hydrodynamique  du  Dorado  devrait  permettre  d’obtenir  une  plus  grande 
precision  avec  le  SIMRMS.  Ce  dernier  est  un  important  outil  de  simulation  pour  la 
planification  et  L  execution  des  operations  de  deminage  a  distance  et  de  leve  des  routes  de 
navigation. 

Plans  pour  Pavenir 

II  est  souhaitable  que  les  algorithmes  ameliores  relatifs  aux  parametres  de  commande  soient 
valides  par  des  essais  en  vraie  grandeur  avec  le  Dorado.  Ces  essais  sont  prevus,  mais  il  n’y  a 
pas  encore  eu  de  fonds  consacres  a  leur  execution. 
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Introduction 

The  Canadian  Remote  Minehunting  Sytsem  (CRMS)  technology  demonstrator  consists  of  a 
remote  drone  (called  Dorado,  8.2  m  long)  towing  a  towfish  (3  m  long)  at  depths  to  200  m. 
The  towfish  carries  a  sonar  which  is  used  for  surveying  and  minehunting  purposes.  The 
CRMS  promises  to  be  a  low  cost,  low  risk  way  of  minehunting  (Figure  1). 

To  help  evaluate  the  system  and  examine  maneuvering  strategies,  DRDC  Atlantic  has 
developed  the  SIMRMS  simulator  [1]  which  provides  a  six  degree-of-freedom  dynamic 
computer  simulation  of  the  two  vehicles,  the  tow  cable  connecting  them,  and  the  control 
systems  they  use.  The  equations  of  motion  in  this  simulation  require  a  description  of  the 
hydrodynamic  forces  on  the  vehicle.  This  description  is  provided  primarily  in  the  form  of 
hydrodynamic  coefficients  (derivatives)  but  some  rational  flow  modelling  can  also  be  used. 

This  report  is  concerned  with  estimating  the  hydrodynamic  forces  for  the  proposed  Build  3 
version  (Fall  2002)  of  Dorado.  The  Build  3  version  is  assumed  to  include  the  new,  enlarged 
fore  and  sternplanes  [2],  but  a  version  of  the  forces  is  also  provided  for  the  original  Dorado 
(Build  1)  fore  and  stemplane  geometry.  The  report  is  not  concerned  with  propulsion  forces 
or  the  hydrodynamics  of  the  tow  cable  or  the  towfish. 

The  hydrodynamic  forces  are  initially  estimated  as  derivatives  by  the  DERIVS  program  [3] 
which  calls  the  DREA  Submarine  Simulation  Program  (DSSP20)  [4]  to  obtain  estimates  of 
the  steady  state  hydrodynamic  forces.  The  ESAM  program  (Estimate  Submarine  Added 
Masses)  [5]  provides  estimates  of  the  acceleration  coefficients  (added  mass  effects). 

While  DSSP20  is  the  basis  for  the  steady  state  derivatives,  its  calculations  require  corrections 
in  two  areas.  First,  it  is  unable  to  predict  the  drag  on  the  unconventional  Dorado  hull  with  its 
many  fittings,  protrusions,  and  holes.  Secondly,  control  deflection  efficiency  in  the  current 
version  (June  2001)  is  incorrectly  modelled.  In  addition  to  these  corrections,  the  Dorado 
snorkel  mast  with  its  self  aligning  fairings  cannot  be  handled  at  all  by  DSSP20.  These 
effects  are  calculated  manually. 

The  mast  effects  are  particularly  challenging  for  two  reasons.  First,  all  the  mast  forces  are  a 
function  of  vehicle  depth  (ie,  the  wetted  length  of  the  mast).  The  snorkel  mast  extends 
through  the  free  surface  and  experiences  considerable  drag,  both  along  its  length  and  right  at 
the  free  surface  (wave  and  ventilation  drag).  So,  for  example,  the  pitch-up  moment  caused 
by  the  mast  varies  with  depth,  something  that  has  been  observed  by  the  operators.  Note  that 
the  wave  drag  modelled  here  is  wave  drag  from  the  mast,  not  from  the  main  vehicle.  This 
latter  effect  is  thought  to  be  small  (the  hull  is  two  to  three  diameters  below  the  free  surface) 
and  it  is  not  modelled  at  the  present  time. 

The  second  reason  the  mast  is  a  challenging  appendage  to  model  is  that,  although  its 
self-aligning  fairings  do  not  generate  lift  in  a  translating  flow,  they  do  generate  lift  when  the 
vehicle  turns.  This  lift  results  from  the  curved  streamlines  interacting  with  the  symmetrical 
fairing.  The  lift  is  reversed  and  amplified  by  the  self-aligning  mechanism,  causing  the 
vehicle  to  roll  to  the  outside  of  the  turn. 
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DERIY S/DSSP20  Steady  State  Predictions 

DSSP20  builds  a  hydrodynamic  model  of  the  flow  around  a  streamlined  vehicle  based  on 
an  input  file  describing  nothing  but  vehicle  geometry.  The  input  file  for  Dorado  is  shown 
in  Appendix  A  and  its  format  is  described  by  Mackay  [4].  This  geometry  is  shown  below 
in  Figure  2: 


Fig.  2a:  Dorado  DSSP20  Geometry  (No  Mast),  Elevation  View 


Fig.  2b:  Dorado  DSSP20  Geometry  (No  Mast),  View  From  Below 
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Fig.  2c:  Dorado  DSSP20  Geometry  (No  Mast),  View  From  Astern 


y 


DSSP20  deals  with  the  vehicle  component  by  component  (hull,  starboard  sternplane,  port 
sternplane,  etc.)  and  each  component  is  shown  as  a  different  color  in  these  figures.  Only 
the  overall  component  geometry  is  important  so,  for  example,  modelling  the 
hemi- spherical  cap  at  the  nose  of  the  hull  is  unnecessary.  Note  the  absence  of  the  mast. 

The  foreplanes  in  Figure  2  deflect  together  to  provide  depth  control.  The  vertical  stabilizer 
above  the  hull  at  the  aft  end  is  fixed  and  houses  the  diesel  engine  exhaust.  The  rudder  is 
the  single  appendage  below  the  tail.  The  stemplanes  provide  both  pitch  and  roll  control. 
The  fore  and  sternplanes  shown  here  double  the  depth,  pitch,  and  roll  control  authority 
relative  to  the  Build  1  design  [2]. 


The  DERIVS  derivs  .  dat  input  file,  described  in  detail  by  Watt  [3],  is  short: 


dorado 

no,  do  not  zero 
8.2  1025.  5.144 
defv  ds :  1  1.0 
def h  do :  1  -1.0 
defv  df :  3  1.0 
def h  dr :  5  1.0 


residual  drag 

2  -1.0 
2-1.0 
4-1.0 


It  gives  the  root  file  name  for  DSSP20  input  and  output  files,  direction  on  how  to  handle  a 
truncated  stem,  the  vehicle  length,  fluid  density,  and  speed  (all  in  SI  units),  and  then  four 
lines  defining  vertical  and  horizontal  plane  control  deflections.  The  ds,  do,  df,  and  dr 
deflections  define  how  pitch  control  (uniform  sternplane  deflection  ds),  roll  control 
(differential  sternplane  deflection  do),  depth  control  (uniform  foreplane  deflection  df),  and 
yaw  control  (rudder  deflection  dr)  are  implemented  by  appendages  1  through  5  in  the 
DSSP20  input  file. 

DERIVS  outputs  a  set  of  hydrodynamic  coefficients  (first  and  second  order  derivatives)  in 
the  body  fixed  axes  shown  in  Figure  2,  except  the  origin  is  shifted  aft,  on  the  hull 
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centerline,  approximately  opposite  the  hull  center  of  buoyancy.  Feldman  [6]  shows  how 
these  derivatives  are  nondimensionalized  and  used  in  the  equations  of  motion. 

In  the  next  section,  the  DERIVS  output  is  corrected  and  enhanced  to  better  model  drag, 
control,  and  the  lift  and  drag  effects  of  the  mast.  The  fully  corrected  steady  state 
hydrodynamic  derivatives  are  presented  in  Appendix  B. 


Correcting  the  Steady  State  Derivatives 

The  DERIVS  generated  derivatives  are  adjusted  to: 

•  correct  vehicle  drag, 

•  correct  appendage  control  derivatives, 

•  correct  foreplane  effectiveness, 

•  include  effect  of  mast  drag  as  a  function  of  depth  and  vehicle  orientation, 

•  include  side  force  generated  on  the  mast  self-aligning  fairings  in  curved  flow  (ie,  when 
the  vehicle  turns)  also  as  a  function  of  depth. 

Mastless  Drag  Correction 

Mastless  drag  is  based  on  full  scale  trials  with  the  DOLPHIN  Mk  1  vehicle  [7].  In 
these  trials,  the  DOLPHIN  vehicle  was  driven  with  an  electric  motor  with  and  without 
the  mast  on.  With  the  vehicle  travelling  at  10  knots,  the  mastless  drag  (hull  plus  all 
appendages  except  the  snorkel  mast)  was  542  lb.  This  gives  a  dimensionless  mastless 
axial  force  coefficient: 


X*  =2 

experiment,  Mk  1,  no  mast 


drag 
pi2  U2 


=  -.003288 


The  DERIVS  estimate  for  the  Mkl  X*  coefficient  is: 


Y*  =  -  001 77 

DERIVS,  Mk  1,  no  mast  ‘ w  J 

The  difference,  A  X*,  is  due  to  the  abnormal  shape  of  the  Mkl  hull,  its  many 
attachments,  external  latch  handles,  holes  and  slots  in  the  hull,  external  weldments, 
etc: 


A  A*  =  -.001558 

This  difference  is  simply  added  to  the  Dorado  X*  coefficient  predicted  by  DERIVS. 

Tailplane  Control  Derivative  Corrections 

DSSP20  (June  2001version)  predicts  tailplane  effectiveness  (the  contribution  of  the 
tailplanes  to  the  transverse  force  when  the  vehicle  is  at  incidence)  correctly  but 
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overpredicts  tailplane  control  derivatives.  It  assigns  an  all-moving  sternplane  the  same 
lift  increment  for  a  control  8v  deflection  as  for  a  vehicle  angle  of  attack  a  deflection. 
In  the  former  case,  only  the  stemplanes  deflect  relative  to  the  flow  while  in  the  latter 
both  hull  and  stemplanes  deflect.  The  former  case  generates  substantially  less  lift  than 
the  latter  case.  And  there  is  a  further  substantial  lift  reduction  required  to  correctly 
model  differential  deflections  (for  roll  control).  For  pitch  control,  equal  lift  on 
simultaneously  deflected  stemplanes  reinforces  itself  across  the  body,  but  the  opposing 
lift  from  differential  deflections  negates  itself  at  the  center. 

Pitts  et  al  [8]  use  slender  body  theory  to  show  that  the  lift  generated  on  a  slender 
cylindrical  body  and  wing  combination  due  to  incidence  a  is: 

Ax  -  *  Av (B)  +  A(WP  Av 

where: 

Av(fi)  +  Avw) -  *  1  +  A 


and  Lw  is  the  lift  on  the  exposed  wing  area  with  body  removed  and  exposed  wings 
joined  at  the  centerline,  KW(B)  is  a  dimensionless  coefficient  giving  the  lift  on  the 
wing  in  the  presence  of  the  body,  and  KB(W)  gives  the  carry-over  lift  on  the  body  in 
the  presence  of  the  wing.  The  X  parameter  is  the  ratio  of  the  hull  diameter  d  to  wing 
tip-to-tip  span  b  (for  the  body  and  wing  combination)  so  that  0  <  X  <  1. 


Pitts  et  al  also  quote  the  results  for  lift  due  to  deflection  8  and  these  can  be  simplified 
to: 

As =  Av(fl)  Av 

so  that: 


T  IS 

iVB(W) 

1  +  K 

W (B) 


Analytical  expressions  derived  from  slender  body  theory  are  available  for  KW(B)  and 
KB{W)  and  are  presented  by  Pitts  et  al.  However,  it  is  approximately  true  that: 


TS 

iVB(W) 

is 

YW (B) 


■  X  +  X  ( 1  —  X)  +  higher  order  terms 


To  first  order,  this  ratio  is  just  X  which  exactly  matches  the  slender  body  theory 
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expression  at  X  =  0  and  1.  The  quadratic  term  in  X  on  the  RHS  is  a  second  order 
(small  perturbation)  correction.  Coupling  this  expression  with  that  for  the 
KW(B)  +  KB(W)  given  above,  it  is  asymtotically  consistent  that: 


Similarly: 


tfW(fl)=l+A,-^A,(l-A.) 

1  1  7i(l -X) 

0+^5  (1+^)2 


As  the  appendage  span  b  goes  to  infinity  ( X  ->  0),  there  is  no  difference  between  8 
and  a  deflections;  the  hull  is  inconsequential.  But  as  the  span  becomes  increasingly 
small  (b  ->  d,  X  ->  1),  a  8  deflection  generates  only  half  the  lift  that  is  generated  when 
the  hull  deflects  as  well. 

Sternplane  Pitch  Control 

Let  C  be  the  correction  that  needs  to  be  applied  to  the  DERIVS  sternplane  pitch 
control  (Ss)  derivatives  (s' for  sternplane,  u  for  uniform  deflection).  Since  these 
derivatives  erroneously  use  Lg  =  La,  they  are  corrected  by  multiplying  by: 

(  l  \ 


That  is: 

We  get: 


1  1  A  (1  -  ^) 


l+X  5(1+  X)2 


C  =C  U 

0  su  o 


correct 


DERIVS 


b  =  2.53708,  d  =  . 96844,  A  =  .3817 
C  =.6990 


where  X  is  calculated  using  the  hull  diameter  d  at  the  root  leading  edge.  Recent 
experiments  suggest  this  as  a  way  to  empirically  correct  for  the  boundary  layer. 

This  correction  is  applied  as  is  to  Zg  and  M. .  The  remaining  Sv  control  derivatives 

S  S 

are  second  order  X  and  M  derivatives.  They  result  primarily  from  induced  drag 
which  increases  as  the  square  of  the  lift  on  the  appendage.  Therefore,  the  X  and  M  Sy 

2 

derivatives  are  corrected  by  multiplying  by  C  . 
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Sternplane  Roll  Control 

There  are  two  further  effects  that  need  to  be  accounted  for  with  differential  deflections. 
First,  the  rolling  moment  from  each  half  of  a  simple  low  aspect  ratio  wing  deflecting 
uniformly  across  its  centerline  (so  that  the  moments  oppose  each  other  and  the  net 
moment  is  zero)  is  twice  that  when  each  half  deflects  differentially  (when  the  moments 
have  the  same  sign;  see  De  Young  [9]).  This  is  because  the  lift  (pressure  distribution) 
for  uniform  deflection  reinforces  itself  across  the  centerline  and  is  actually  greatest 
there,  whereas  the  lift  for  differential  deflections  cancels  at  the  centerline. 

Second,  with  uniform  sternplane  deflections,  the  velocity  field  that  is  generated  is 
symmetrical  about  the  vertical  centerplane  through  the  hull.  So  there  is  no 
intereference  with  the  vertical  tailplane  appendages.  For  differential  sternplane 
deflections,  however,  crossflow  is  present  at  the  centerplane  (resulting  from  sternplane 
tip  and  strong  hull  bound  vorticity)  which  generates  a  moment  on  the  vertical 
appendages  counteracting  the  moment  developed  by  the  sternplanes  themselves 
(Adams  and  Dugan  [10]). 

So  there  are  three  good  reasons  why  roll  control  is  best  implemented  with  ailerons 
located  well  outboard  of  the  hull  centerline: 

•  differential  lift  (opposing  pressure  distributions)  is  not  compromised  by  close 
proximity  of  the  lifting  components, 

•  crossflow  on  the  centerplane  is  minimized,  thereby  minimizing  roll  reduction  from 
interaction  of  the  crossflow  with  vertically  oriented  appendages  on  the  centerplane, 
and 

•  rolling  moment  is  maximized  because  of  the  larger  moment  arm. 

For  these  reasons,  it  makes  good  sense  to  implement  roll  control  on  the  Dorado  mast 
rather  than  its  sternplanes.  Doing  so  provides  the  added  benefit  of  allowing  the 
sternplanes  to  be  devoted  entirely  to  pitch  control.  However,  roll  control  on  the  mast 
requires  an  additional  system  and  this  prevents  its  consideration  at  this  time. 

De  Young  [9]  uses  slender  body  theory  to  analyse  the  roll  control  generated  by 
differentially  deflected  components  on  an  isolated  wing.  If  the  deflecting  components 
are  symmetrical  about  the  centerline  and  are  the  all-movable  ends  of  the  wing,  and  A 
is  the  ratio  of  the  nondeflecting  inner  portion  of  the  span  length  to  the  overall  span 
length,  then  the  dimensionless  rolling  moment  generated  by  differentially  deflecting 
each  wing  component  an  angle  8o  is: 


3 


W(A,A)  =  2 


moment 


p  U2Sb 
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■a-*-) 


(3/2) 


2 

where  W(A.,  A)  is  a  ‘wing’  alone  function,  S  =  be  is  the  wing  surface  area,  A  =  b"/S  = 
blc  is  the  aspect  ratio  of  the  wing,  and  c  is  the  mean  chord  length. 

When  the  nondeflecting  inner  portion  is  replaced  with  a  circular  body,  so  the 
deflecting  ends  of  the  wing  are  now  sternplanes,  Adams  and  Dugan  [10]  use  slender 
body  theory  to  show  that  the  dimensionless  moment  WB(X  A)  (for  ‘wing  and  body’) 
is  a  complicated  numerical  integration  of  an  elliptic  function.  Their  plotted  results  can 
be  fit  with: 


WB(A.,  A)  = 


1  ,  (3/2)  7  in  7  7 

7  (1  -  X2)  +  .720  X2  (1  -  X)10  +  -819  X2  (1  -  Xy 

6 


.520  X 3  (1  -X)2-  1.306  X12  (1  -X)2 


Adams  and  Dugan  also  include  the  effect  of  the  vertical  fin  (having  the  same  planform 
as  the  sternplanes)  on  the  moment  generated  by  differentially  deflecting  sternplanes: 


WBF(X,  A)  =  (.127  (1  -  X2)3  /  "  ’  +  2.043  X2  (1  -  X)2  -  1.979  X2  (1  -  X)4 
-  2.543  X3  (1  -  X)2  -  .661  X13  (1  -  X)2)A 

These  three  functions  are  shown  in  the  next  figure. 
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Fig.  3:  Roll  Control  with  an  Isolated  Wing  (W), 
a  Wing  and  Body  (WB),  and  a  Wing,  Body,  and  Fin  (WBF) 


Adding  a  body  to  the  wing  apparently  has  two  major  effects.  First,  it  helps  to  isolate, 
and  therefore  maintain,  the  opposing  loads.  This  effect  tends  to  increase  the  moment 
and  it  dominates  for  small  body  diameters  when  the  opposing  loads  are  close  to  each 
other.  Second,  it  reduces  the  amount  of  lift  supported  by  the  adjacent  portion  of  the 
inner,  nondeflecting  part  of  the  wing.  This  effect  is  only  seen  for  large  nondeflecting 
inner  span  lengths,  when  X  approaches  1  and  the  opposing  loads  are  at  the  ends  of  the 
wings  well  away  from  one  another. 

Adding  a  vertical  fin  reduces  the  moment  since  the  fin  resists  the  circulation  field  the 
deflecting  stemplanes  generate. 

Adams  and  Dugan  note  that  if  the  vertical  fin  also  sustains  a  differential  deflection  5o, 

then  the  net  moment  generated  is  2  WBF(X,  A).  In  linearized  slender  body  theory,  the 
principle  of  superposition  applies. 

Dorado’s  vertical  appendages  are  smaller  than  its  stemplanes,  so  the  vertical 
appendages  counteract  the  sternplane  induced  roll  less  than  if,  as  the  WBF  function 
assumes,  they  had  the  same  planforms  as  the  stemplanes.  However,  this  decrease  in 
the  reduction  is  assumed  to  be  negligible  since  the  missing  span  length  is  well  away 
from  the  body  and,  therefore,  also  the  sternplane  induced  circulation  impinging  on  the 
vertical  fins. 

DSSP20  assumes  that  roll  control  is  simply  the  lift  that  would  be  generated  on  a  single 
appendage  with  the  vehicle  at  incidence  times  4/ (3k)  of  the  half  span  length  (the 
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spanwise  center  of  pressure  for  half  the  lift,  assuming  an  elliptic  lift  distribution)  times 
2,  since  there  are  two  sternplanes.  This  calculation  accounts  for  the  effect  of  the 
boundary  layer  and  afterbody  contraction  on  tailplane  efficiency.  The  correction  that 
needs  to  be  applied  to  the  DSSP20  prediction,  therefore,  is  the  ratio  of  the  slender  body 
theory  prediction  of  actual  roll  control  moment  (WBF(3,,A))  to  that  from  incidence  lift 

acting  at  the  4/(371)  spanwise  location.  This  latter  moment  is  KW(B)  +  K[t>(  w  ( =  (1+3,)' 
times  the  slender  body  theory  prediction  of  dimensionless  lift  on  half  the  exposed  wing 
(7t  a/4  where  a  is  the  aspect  ratio  of  the  exposed  sternplanes  joined  together),  times 
S  /  S.  to  convert  the  nondimensionalization  of  this  latter  lift  from  an  area  based  on  the 

a  A 

exposed  sternplanes  to  one  based  on  the  full  span  planes,  times  the  dimensionless 
moment  arm  4/(37t)/2,  times  2  because  there  are  two  sternplanes: 

1  (l+X)2aSa 

3 


So,  the  correction  Csd  (cl  for  differential)  to  be  applied  to  DERIVS  8  derivatives 
depends  only  on  X : 


2  x  moment  from  incidence  lift  on  one  sternplane  = 

1  9  2 

=  ~A(l-X2) 

since: 

aSo  _(b-df 
ASa=  b* 


c 


sd 


WBF(X,A) 
A(1  -X2) 


=  .4845 

2 

Csd  is  applied  as  is  to  derivatives  linear  in  8o  and  Csd  multiplies  those  second  order 
in  8  2. 

O 


Rudder  Control 

The  single  rudder  appendage  deflects  opposite  the  stationary  vertical  stabilizer  and  in 
the  presence  of  the  sternplanes.  From  a  slender  body  theory  perspective,  this  is  a 
superposition  of  two  effects:  the  rudder  and  vertical  stabilizer  deflecting  uniformly 
8/2  degrees  and  the  two  planes  deflecting  differentially  8/2  degrees. 

The  uniform  8/2  deflection  determines  the  side  force  and  yawing  moment  response. 
This  will  be  half  the  response  of  uniformly  deflecting  symmetrical  rudder  appendages 
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deflecting  8.  degrees.  DSSP20  tries  to  do  this  but  errs  in  predicting  the  uniform  load 
in  the  same  way  as  for  uniform  sternplane  deflection.  The  correction  is,  as  before,  to 
multiply  the  DERIVS  Fg  and  N~  derivatives  by  a  correction  Cru  where: 

r  r 


c 


fi  \ 


L 

V  “7 


1  1  X(l  -?i) 

~  l +X  5  (1  +  X)2 

=  .6281 

where: 

b  =  1.646,  d  =  .89,  *,=  .5407 

with  d  again  the  local  hull  diameter  opposite  the  leading  edge  root  location. 

The  differential  8/2  deflection  primarily  determines  the  rolling  moment  derivative 
K& .  The  derivative  is  half  the  prediction  for  symmetrical  rudder  appendages 

r 

differentially  deflecting  8  degrees.  Again,  DSSP20  tries  to  do  this  but  errs  as  above 
for  the  differential  sternplane  deflections.  The  correction  is  to  multiply  the  derivative 
by  Crd  where: 


C 


rd 


=  3 


WBF(/i,  A) 
A(  1  -X2)' 


=  .5452 


Actually,  the  Ks  derivative  also  has  a  contribution  from  the  Fg  derivative  since  the 

r  r 

Dorado  stem  droops.  That  is,  the  center  of  the  imaginary  symmetrical  pair  of  rudders 
uniformly  deflecting  8/2  degrees  is  0.077  m  below  the  centerline  of  the  hull  (at  the 
rudder  leading  edge)  about  which  the  rolling  moment  is  calculated.  So  the  DERIVS 
K5  derivative  is  actually  calculated  from: 

r 


(2  b  \ 
~  ~  +  .077 
[3% 


so  that  0.011 1 (0.011 +2b/3n)  of  the  DERIVS  7fg  derivative  is  attributable  to  uniform 

r 
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loading  and  therefore  requires  the  Cru  correction,  and  the  remainder  is  attributable  to 
differential  loading  and  requires  the  C  ,  correction.  Therefore,  the  net  correction  to 
is: 


r 


CrK=  -1806  Cn  +  .8194  CM 
=  .5602 


which  is  not  very  different  from  the  previous  calculation. 

The  DERIVS  Tg  and  Vg  derivatives  are  corrected  by  multiplying  by  Cyu.  The  V. 


r 


r 


r 


derivative  is  corrected  by  multiplying  by  CrK.  The  X  and  M  second  order  8. 

derivatives  are  inherently  nonlinear  and  no  linear  theory  or  superposition  argument  can 
explain  how  to  correct  them.  A  second  order  correction  that  is  some  kind  of  average 
of  the  Cru  and  CrK  correction  terms  seems  appropriate,  however,  and  so  we  use  the 

simple  product  Cm  CrK  to  correct  the  second  order  derivatives. 

Foreplane  Control  Derivative  Corrections 

The  DERIVS  estimated  foreplane  control  derivatives  are  known  to  be  underestimated. 
In  this  section,  we  recalculate  foreplane  effectiveness  to  see  where  the  error  occurs. 

The  methods  of  Pitts  et  al  [8]  again  form  the  basis  of  the  analysis. 

The  lift  generated  by  the  foreplanes  due  to  an  angle  of  attack  a  is  modelled  via: 


as  described  in  the  previous  section.  The  lift  L w  is  the  lift  on  the  exposed  appendages 

with  the  body  removed  and  the  appendages  joined  at  the  centerline  (the  ‘wing  alone’). 
This  lift  is  estimated  using  Whicker  and  Fehlner’s  [11]  expression  for  the  lift  on  a  low 
aspect  ratio  isolated  wing: 


d 


+  4 


where: 


L 


p  U2S 

r  a 


and  £1  is  the  sweepback  of  the  quarter-chord  line.  This  gives  the  foreplane 
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contribution  to  the  normal  force  derivative  Z  as: 

W 


A  Z  =— 2 

W 


p  f/2/2 


But: 

so: 


Since: 


(1  +XfnaS 


-1.8 


r 


a 


1.8+ cos (£2),.,  4 

V  V  COS(Q) 


+  4 


aSa  =  b2(l  -X)2 


A  Z  =-1.8 


(1-X2)  nb2 


1.8  +  cos(£2)  < 


cos(£2) 


'  +  4 


Z?  =  2.5370  m 
a  =2.484 
^  =  .4605 
£2  =  15.27  degrees 

the  foreplane  contribution  to  Z  is: 

r  w 


A  Z  =-.06697 

W 

The  DERIVS/DSSP20  prediction  is  17%  larger  than  this,  which  seems  large  but  is 
probably  within  the  error  of  these  estimates. 

From  the  previous  section  we  know  that: 

h_  1  1  7i(l -X) 

4  _  1  +  ^  5  (1+4* 


so  that: 


=  .6614 
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1  1A(1-A) 


f 


1  +x 


\ 


5  (1+A)2 


A  Z 


=  -.04430 

And  this  is  20%  larger  than  the  DERIVS/DSSP20  prediction  of  Zg  .  That  is,  the 

effective  DSSP20  prediction  for  Lg/  La  is  0.47  which  is  too  low  since  we  know  this 
ratio  achieves  a  lower  limit  of  0.5  as  A,  ->  1.  A  much  larger  underprediction  of  Zg 
occurs  with  the  original  Dorado  foreplanes,  as  shown  in  Appendix  D. 

Two  sets  of  corrections  are  made.  First,  Zw,  Mw,  Z  ,  and  Mq  are  adjusted  by  adding 
the  following  DZw,  DMw,  Dz  ,  and  DMq  values  to  account  for  the  discrepancy  in  Zu>: 


D7  =  -.06697  -  (A  Z  ) 

Zw  v  W/ 


DSSP20 


=  .01120 


D7  ( x  -x„„) 

Zw  v  qc  CB 7 

D..  =- - 

Mw  1 

=  -.002322 


D7  =  D. . 

Zq  Mw 

=  -.002322 


D7  (. x  -  x„D) 

Zq  v  qc  CB 7 

- 7 - 

=  .0004815 


where  x  -  xCB  is  the  distance  from  the  center  of  buyancy  to  the  nominal  quarter 
chord  location  of  the  foreplanes. 

No  change  is  made  to  drag  or  second  order  derivatives  as  a  result  of  these  corrections. 
The  corrections  are  small  and  may  be  explainable  by  the  use  of  different  algorithms. 
Nevertheless,  they  are  made  in  order  that  the  derivatives  are  consistent  with  the 
algorithms  used  to  design  the  foreplanes  [2]. 

The  second  correction  is  to  Zg  which  is  achieved  using  the  factor: 
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1 


c  =-.04430 - 

8/ 

DSSP20 

=  1.203 

The  linear  derivatives  are  corrected  by  multiplying  by  .  The  second  order  8; 

2 

derivatives  are  corrected  by  multiplying  by  Cf  . 

Mast  Drag 

The  mast  forces  are  a  function  of  the  length  of  the  submerged  mast,  sm  .  This  varies 
with  the  depth  of  the  vehicle  and  its  pitch  and  roll  angles: 


Fig.  4:  Mast  Coordinates  and  Vehicle  Pitch  and  Roll  Euler  Angles 


View  Along 
Body  x  axis 


_  %  ~  Xm  Sill(9  ) 

S"'  cos(0)cos((]))  r,n 

where  zt)  is  an  inertial  coordinate  giving  the  depth  of  the  origin  of  the  vehicle  body 
fixed  axes,  x  is  the  axial  coordinate  of  the  mast  in  body  fixed  axes,  0  and  <f»  are  the 
vehicle  pitch  and  roll  Euler  angles,  and  r  is  the  hull  radius  at  the  mast: 
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body  axes  origin  =  3.7  m  aft  of  nose 
x  = .4234  m 

m 

r  =  .5842  m 

m 

Mast  drag  is  based  on  full  scale  trials  with  the  DOLPHIN  Mk  I  vehicle  [7].  The 
DOLPHIN  and  Dorado  self-aligning  fairings  are  geometrically  similar,  although 
submerged  mast  length  and  other  vehicle  properties  are  not.  The  Dorado  mast  chord  is 
largest: 


Dorado  mast  chord  length 

- =  1.209 

DOLPHIN  Mk  I  mast  chord  length 

and  the  square  of  this  ratio  is  used  to  scale  experimentally  measured  DOLPHIN  forces 
to  Dorado  forces  (except  for  the  section  profile  drag  which  scales  in  proportion  to  the 
ratio). 

In  the  trials,  DOLPHIN  was  powered  by  an  electric  motor  with  and  without  the  mast 
attached,  allowing  the  mast  drag  to  be  isolated.  With  the  vehicle  travelling  at  10  knots, 
the  mast  drag  was  353  lb.  The  component  contributions  to  this  drag  are  estimated  as 
[7]: 


•  Ventilation  drag  (a  vent  a  foot  or  so  deep  opens  up  at  the  surface  behind  the  fairings) 
=  120  lb. 

•  Spray  drag  =  13  lb. 

•  Wave  drag  =  130  lb. 

These  three  drag  components  act  at  the  surface  at  (xm,  z)  in  body  axis  coordinates. 

They  are  totaled,  scaled  to  Dorado,  and  nondimensionalized  (using  the  square  of 
Dorado’s  hull  length)  in  d  ,  the  surface  component  of  mast  drag: 

ds  =  .001871 

zs  =  -rm~sm 


•  Junction  drag  =  13  lb. 

This  is  scaled  and  nondimensionalized  in  d,,  which  acts  at  the  base  of  the  mast  at 

K,  zb ): 


db  =  .00009250 
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•  Profile  drag  =  77  lb. 

This  is  converted  to  a  dimensionless  section  profile  drag  d  (drag  per  meter)  using 
the  fact  that  the  nominal  DOLPHIN  submerged  mast  length  is  9.5  feet. 

d  =.0001565 

pm 


The  Dorado  profile  drag  is  reconstituted  by  integrating  d  along  the  mast  in  the 
presence  of  a  local  velocity  field. 

Because  the  fairings  self  align,  the  drag  acts  locally  in  the  direction  of  the  local  flow 
(assuming  there  are  many  separate  fairings  in  the  spanwise  direction,  which  is  not 
strictly  correct).  Mast  drag  is  generated  by  5  of  the  6  possible  vehicle  velocities: 
u,  v,  p,  q,  r  (see  Watt  [5]).  The  first  two  are  the  axial  and  lateral  translational 
velocities;  the  third  translational  velocity  w  (co-axial  with  the  mast),  is  assumed  to 

not  influence  drag.  The  roll,  pitch,  and  yaw  rotational  velocities  p,  q,  r  do  contribute 
to  the  drag.  The  influence  of  all  5  motions  on  the  drag  is  a  function  of  their  integrated 
effect  along  the  mast  in  a  local  flow  influenced  by  the  hull. 

The  axial  velocity  of  the  mast  um  is  a  function  of  the  vehicle  axial  velocity  u,  its 
pitching  moment  q,  and  the  position  along  the  mast  z  : 


um  =  u+qz 


There  is  no  interference  with  axial  force,  assuming  the  hull  is  an  infinitely  long 
cylinder. 


The  lateral  velocity  of  the  mast  v  depends  on  the  vehicle  lateral  velocity  v,  its  yaw 
angular  velocity  r,  its  roll  angular  velocity  p,  the  axial  location  of  the  mast  x  ,  the 
local  hull  radius  at  the  mast  rm,  and  the  location  along  the  mast  z  : 


v  =  (v  +  r  x  ) 

m  m/ 


1  + 


2  A 


2 

*  7 


-pz 


The  lateral  velocity  of  the  hull  opposite  the  mast,  v  +  rxm,  generates  an  interference 
velocity  along  the  mast  which  is  modelled  using  the  velocity  field  from  2D  potential 
flow  over  a  circular  cylinder.  The  lateral  velocity  generated  by  roll,  -p  z„  generates 
no  interference  velocity  about  the  axisymmetric  hull. 

Since  crossflow  velocity  was  not  present  in  the  trials  from  which  the  mast  drag 
components  were  extracted,  the  mast  drag  coefficients  are  correct  as  shown  above. 

With  Dorado  moving  with  6  degrees-of-freedom,  the  local  dimensionless  drag  force  on 
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2  2  2  2  2  2  2 

the  mast  is  d.  (u  +v  )/U~  where  U  =  u ~  +  v  +w  .  The  axial  and  lateral 

i  v  m  m  7 


components  of  this  force  are  obtained  by  multiplying  by  the  cosine  u 


/  /  2  2 

/  A/  w...  +  v_ 


m  m 


and  sine  v 


A fu 


m  m  m 


2  2 

+  v  of  the  crossflow  angle: 


/oca/  axial  drag  force  component - 


2  2 

d  u  \l  u  +  v 

i  m  V  fw 


u1 


local  lateral  drag  force  component - 


2  2 

d  v  \l  u  +v 

i  m  V  m  m 


uz 


Of  course,  these  forces  act  in  a  direction  opposite  to  the  w  and  v  mast  velocities. 

5  rr  mm 

It  is  convient  to  define  axial  Fax(z)  and  lateral  F { t(z)  functions  from  these  force 
components  for  calculating  the  forces: 


F  ( z)  =  u 

riY'-''  n 


2  2 

+  V 


(u  +  q  z)/\  I  (u  +  qz)  + 


f 

f  2  A 

r 

m 

\ 

(v  +  rxj 

l  +  — 

-pz 

V 

l  z  ) 

) 

V2  2 

Um  +V,n 


f 

f  2  A 

r 

m 

\ 

L  * 

f 

f  2  ^ 

r 

m 

y 

(v  +  rx  ) 

v  my 

1  +~ 

-pz 

/\  (u  +  qz)  + 

(v  +  rx  ) 

v  m7 

1  +~ 

-pz 

V 

l  z  ) 

V 

V 

l  Z  ) 

y 

So,  when  calculating  the  forces  at  the  surface,  for  example,  the  axial  and  lateral  force 
components  are  proportional  to  ds  Fafzs )  and  ds  Flat(zs).  For  the  base,  the  b 
subscript  is  used.  However,  to  get  the  profile  drag,  d  F.(z )  needs  to  be  integrated 
along  the  length  of  the  mast.  And  to  get  the  moment  generated  by  the  profile  drag,  we 
need  to  multiply  the  integrand  by  z.  This  gives  rise  to  the  four  integrals  calculated 
below. 
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Unfortunately,  integrals  of  Fax  and  F/at  cannot  be  solved  analytically.  And  it  is 

unwieldy  to  have  to  do  a  numerical  integration  every  time  step  of  a  simulation.  So  a 
conventional  expansion  of  these  functions  is  made  before  they  are  integrated.  The 
expanded  and  truncated  functions  are  designated  f  (z)  and  flat(z).  The  expansion 

assumes  the  axial  velocity  u  is  much  larger  than  any  other  velocity  component,  and 
third  and  higher  order  terms  are  truncated.  Thus,  the  mast  profile  drag  contribution  to 
the  forces  is  second  order  accurate,  which  is  consistent  with  the  other  derivatives 
calculated  by  DERIVS. 

The  expansion  for  the  axial  force  function  F  (z)  and  its  associated  integrals  are: 


(  4  \ 


fjz)  =  u~+2qzu+: 


1 


f 

f  2  A 

r 

m 

\ 

(v  +  rx  ) 

v  m7 

1  +~ 

-pz 

V 

l  ^  ) 

) 

1  = 

ax 


f.,lz)  dz 


2  1  2  2  2 

■u  s  -  ci  u  s  ( s  +2r  )+  a  s  (s  +  3  r  s  +3  r  ) 

m  1  m  v  m  m7  Q  1  m  x  m  mm  m  ' 

1  3  3  2  2  2 

^  s  (3  s  +12  r  +24  r  s  +16  r  s  )(v  +  rx  ) 

in  '  m  m  m  m  m  m  7  nr 

3 

(r  +5  ) 

v  m  m/ 


+  ■ 


+  ^  p  (v  +  rx  ) 

O  1  v  m/ 


r 

f r  +  s  ^ 

s  +2  r  s  +2  r  In 

m  mm  m 

m  m 

V 

\  m  JJ 

1  2  2  2 

+  ^  p  s  (s  +3  r  s  +3  r  ) 

fa*  m  x  m  mm  m  ' 


zfjz)  dz 


I  2  2  2  2 

~2U  S,n  (Sm  +  2  Ym)  +  3  C!  U  S,n  (Sm  +  3  Tm  Sm  +  3  r,n  ) 

1  2  2  2  2 

-  a  s  (s  +  2  r  )  (2  r  +2  r  s  +  s  )-r  (v  +  rx  )  In 

A  1  m  v  m  m/  v  m  mm  m  7  m  v  m7 


7r  +s  3 

m  m 


V  m 
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2  2  2 

1  s  (v  +  rx)  (s+2r)  (2  r  +2  r  s  +  s  ) 

i  m  v  m7  v  m  m7  v  m  mm  m  7 

~  4  2 

4  (r  +s  ) 

v  m  m 7 

1  2  2  1  ,  2  2 

p  s  (s  +  3  r  s  +  6  r  ) (v  +  rx  ) - ~ p"  s  (5  +  2  r  ) (2  r  +2  r  s  +s  ) 

'l  1  m  y  m  m  m  m  7  v  m7  0  1  m  v  m  m7  x  m  mm  m  7 

Similarly,  for  the  lateral  force  function  FAz ): 


2  2 

v  \l  u  +  v  =v  u  +  v  qz+  0 

m  V  m  m  m  m  1 


(  3  '\ 
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/J^): 
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f  2  'N 
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\ 

f 

f  2  ^ 

r 
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\ 

(v  +  rx  ) 

v  m7 

1  +— r 

-pz 

u  + 

(v  +  rx  ) 

v  m7 

1  +_r 

-pz 

V 

l  z  ) 

V 

l  z  ) 

y 

/ 


1  It  o  mj  h  +2r  )(v  +  rx  ) 

1  1  2  2  m  v  m  m7  x  m7 

=  ~  (s  +2  r  )u  p  s  (s  +3  r  s  +3  r  )q  p  s  + 

9  v  m  m7  *  m  'l  v  m  mm  m  7  1  r  m 


r  +s 

m  m 


-  (v  +  rx  ) 

2  v  m7 


r 

(r  +  s  Y| 

s  +2  r  s  +2  r  In 

m  m  m  m 

m  m 

V 

K  m  JJ 
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Z, 

~  b 

zflat(z)dz 

Z 
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1  2 

-  „  up  s  (s  +3r  i  +3 

O  l  mA  m  m  m 


)  +  A  p  q  s 

7  A  ^  1  n 


(s  +2  r  )  (2  r  +2  r  s  +  s  ) 

\  \  m  mi  m  m  ' 


1 


u  (v  +  rx  ) 

\  m/ 


s  +2  r  s  +2  r  In 

m  mm  m 


r  +s 

m  m 


1  2  2 

+  „  s  q  (s  +  3  r  s  +  6  r  )  (v  +  rx  ) 

2  m  1  v  m  m  m  m  7  v  m7 


\\ 

)) 


Note  that  force  derivatives  are  easily  obtained  from  these  integrated  functions  because 
the  terms  are  all  linear  in  products  of  two  velocities.  If  derivatives  for  the  mast  surface 
and  base  forces  are  also  required,  then  those  forces  should  be  calculated  using  the  f 
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and  flat  functions  instead  of  the  Fax  and  Fh  functions.  For  example,  the  mast  drag 
contribution  Xd  to  dimensionless  axial  force  X  is: 


d  u  \  u  +  v 

s  m  V  m  m 


uz 


d,  u  \  u  +v 
b  m  V  m  m 


uz 


1  2 

2 

d  u  \ 

/ 11 

+  V 

pm  m  \ 

j  m 

m 

u1 


dz 


dJaA)  dbfa^b)  dP’"L 

u 2  u 2  u 2 


X,  U1  X,  uql  X,  q2l 2 


U1 


'  +  ■ 


uq 


uz 


■  +  ■ 


m 


Uz 


'  +  .  .  . 


which  allows  the  dimensionless  derivatives  Xd  ,  Xd  ,  etc.  used  by  Feldman  [6]  to  be 

uu  uq 

obtained. 


A  second  order  representation  is  usually  adequate  for  describing  the  hydrodynamics  of 
a  vehicle  for  simulation  purposes  (excluding  extreme  maneuvers).  This  is  thought  to 
be  particularly  true  for  Dorado  since  its  main  purpose  is  to  ensure  the  stability  of  the 
sonar  it  is  towing.  We  check  the  accuracy  of  the  fax  and  f,  truncated  expansions  by 
comparing  them  with  their  accurate  F  and  Ft  parents  in  a  worst  case  situation.  For 
the  mast,  the  expansions  are  least  accurate  where  the  q  z  and  vm  velocities  are  largest. 
This  is  at  the  surface  where  z  =  z  .  During  a  typical  DOLPHIN  10  knot  maneuver  (u  = 

5  m/s),  full  scale  trials  results  show  that  p,q,r  maximum  magnitudes  may  be  as  large 
as  0.3,0.17,0.5  rad/s  in  some  situations.  Combining  these  magnitudes  in  the  worst 
possible  way  in  the  Fax,  Fjat,  fax,  flat  functions  and  assuming  zs  =  -3.5  m,  we  find: 


/  (-3.5) 

J  axx  7 

F  (-3.5) 


1.001 


/J-3-5) 

^J-3-5) 


.9610 


so  that  the  truncated  terms  in  the  /  and  flat  functions  are  not  significant  for  Dorado. 
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Dimensionless  Mast  Drag  Forces 

The  dimensionless  mast  drag  forces  are: 


dsFJ&  db  FA)  dp>»L 

u2  "  u2  u 2 

dsFla&)  db  dP™Il«t 

u 2  u 2  u 2 

z/=o 


dsZsFJ Z,) 

dbZbFJZb) 

d  Iz, , 

pm  lat 

2  "t~ 
U2l 

2  "t~ 
U2l 

ir, 

ds  Zs  F ax'Zs ) 

dbzbFJzb ) 

d  Iz 

pm  ax 

U2l 

U2l 

U2l 

dsXmFJZs) 

dbXmFJZb ) 

d  x  1 

pm  m  lat 

U2l 

U2l 

U2l 

The  X(/  and  Md  expressions  here  give  the  Dorado  mast  drag  at  10  knots  as  481  lb 

(with  z0  =  3  m),  up  from  353  lb  for  the  DOLPHIN  Mkl  vehicle  (zi}  =  2.9  m)  at  the 

same  speed.  The  center  of  action  of  this  drag  is  2.7  m  above  the  hull  centerline,  or  a 
third  of  a  meter  below  the  surface.  This  reflects  the  fact  that  ventilation,  spray,  and 
wave  drag  account  for  about  3/4  of  the  mast  drag  and  they  are  assumed  to  act  right  at 
the  surface. 

An  examination  of  the  other  formulae  shows  that  the  pitching  and  rolling  moments 
vary  at  least  lineary  with  the  mast  wetted  length  s  ,  and  at  least  quadratically  for 
angular  velocity  motions. 

Mast  Lift 

In  translational  flow,  it  is  assumed  that  the  self-aligning  fairings  self-align,  generating 
no  lift  forces  (forces  perpendicular  to  the  local  flow).  In  rotational  flow,  the  fairings 
still  self-align  (ie,  no  moment  is  sustained  at  the  pivot  point,  neglecting  friction),  but 
lift  is  still  generated  because  of  the  effective  ‘camber’  of  the  symmetrical  fairing  in  the 
curved  streamlines.  This  effect  has  been  analyzed  by  Watt  [13]  who  uses  thin  airfoil 
theory  to  show: 

•  The  lift  from  the  curved  streamlines  acts  at  the  midchord  point  x  and  is  directed 
towards  the  center  of  curvature. 

•  The  2D  curvature  lift  coefficient  (based  on  chord  length)  is  C,  =  n  c/2R,  where  c  is 
the  fairing  chord  length  and  R  is  the  streamline  radius  of  curvature. 
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•  These  thin  airfoil  theory  predictions  agree  with  Reynolds  Averaged  Navier-Stokes 
(RANS)  predictions  to  within  8%. 

•  The  curvature  lift  deflects  the  fairing  midchord  towards  the  center  of  curvature, 
creating  an  angle  of  attack  a  which  generates  a  lift  at  the  quarter  chord  (the 
‘aerodynamic  center’  xa )  directed  away  from  the  center  of  curvature  (see  figure  5). 
Equilibrium  is  achieved  when  the  moment  about  the  mast  pivot  point  x  (which 
coincides  with  is  zero. 


Fig.  5:  Streamline  Curvature  Induced  Side  Force  on 
the  Dorado  Self-Aligning  Fairing 


Thus,  the  net  force  on  the  self-aligning  fairing  is  several  times  larger  than  the  original 
curvature  lift  and  in  the  opposite  direction.  Note  that  the  amplification  shown  here  is 
very  sensitive  to  the  correct  determination  of  the  aerodynamic  center  since  the  lift  is 
inversely  proportional  to  (£^  -  qj  which  is  a  small  quantity. 

A  possible  problem  with  self-aligning  fairings  is  that  they  are  subject  to  friction  at  the 
pivot  point.  These  fairings  actually  pivot  about  the  hollow  snorkel  which  has  a  radius 
about  equal  to  t  -  E  .So  the  friction  acts  with  the  same  moment  arm  as  the 

incidence  lift  Cl  .  In  the  past,  it  has  been  assumed  that  friction  would  act  to  oppose 
deflection  by  Cjc  thereby  reducing  the  magnitude  of  C,  .  Comparisons  with  trials 
results  [10]  and  thin  airfoil  theory  predictions  of  Clc ,  qr ,  and  l^ac  suggested  that  an 
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empirical  correction  factor  of  0.8  accounted  for  this  effect.  However,  the  RANS 
predictions  of  these  quantities  [13]  suggest  =  0.46  instead  of  0.5  and  c,ac  =  0.254 
instead  of  0.25  (thin  airfoil  theory  and  RANS  predictions  of  C,  agree  well). 
Although  these  are  only  small  differences,  when  used  in  the  expression  for  C;  below, 
they  eliminate  the  need  for  the  empirical  correction  factor.  Thus,  the  RANS 
predictions  for  these  quantities  are  used  below.  One  assumes,  then,  that  friction  is  not 
important. 

Therefore,  the  net  side  force,  quoted  as  a  2D  lift  coefficient  nondimensionalized  by 
chord  length,  is: 


i 

c'“2  fi(LA„> 

This  is  converted  to  standard  dimensionless  form  (per  meter  of  submerged  mast 
length)  in: 


/,  =  ■ 


Cj  c 


where  f0  is  dimensionless: 


.4 

R 


/o 


C,cR 


jkc2£c-^J 

2 

RANS  calculations  [13]  give  the  following  ^  and  l^ac  values,  allowing  fQ  to  be 
evaluated: 

^=.46  c,  £,ac  =  .254  c,  ^=.1776  c,  c  =  .6223,  /  =  8.2 


/0  =  . 02439 

The  radius  of  streamline  curvature  at  the  mast  is  the  velocity  component  in  the 
horizontal  plane  divided  by  the  angular  rotation  in  the  same  plane: 
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R  = 

r 

where  R ,  and  therefore  /  ,  have  the  sign  of  the  angular  velocity  r. 

2  2  2 

As  with  mast  drag,  we  convert  /  into  a  local  lift  by  multiplying  by  (w  +  v  )  /  U  , 

the  local  dynamic  pressure  contributing  to  the  lift  nondimensionalized  by  the  global 
dynamic  pressure.  This  is  then  multiplied  by  the  cosine  and  sine  of  the  local 
horizontal  plane  flow  angle  to  give  the  lateral  and  axial  force  components: 


/  u 

J  c  m 

local  lateral  lift  force  component  =  — 


local  axial  lift  force  component  = 

JorV,n 

~  u 2 

This  analysis  ignores  the  second  order  induced  drag  from  /  ,  which  is  proportional  to 

2 

the  square  of  the  lift  divided  by  the  aspect  ratio:  0((fQ  r)  /  a),  where  a  is  the  aspect 

ratio.  Although  the  Dorado  mast  has  a  submerged  aspect  ratio  only  of  about  5,  it  has 
the  hull  on  one  end  and  the  free  surface  on  the  other  and  so  effectively  has  a  larger 
aspect  ratio.  This,  together  with  the  fact  that  /0  «  1,  suggests  that  neglecting  the 
induced  drag  term  is  justified. 

The  dimensionless  lift  forces  can  now  be  calculated. 
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i  fns  (s  +  2  r  )  (2  v  r  +  (r  +s  )p  r +2  r2  x  ) 

±  J  l)  m  x  m  my  v  v  m  ra7^  ra7 

~2  f/2(r  +5  ) 

v  m  m/ 

Xt  v  rl  Xj  p  r  l2  Xj  r2 12 

vr  pr  rr 

u2  +  u 2  +  u2 

where  this  last  form  can  be  used  to  identify  coefficients,  if  required.  Since  many  of 
these  coefficients  are  not  now  required  in  a  simulation,  so  that  new  programming 
would  be  required  to  accomodate  them,  it  is  easier  just  to  evaluate  the  entire 
dimensionless  X.  force,  the  second  line  above,  in  a  subroutine  that  calculates  all  the 
mast  effects  (Appendix  C). 

Similarly: 
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The  rotation  lift  generated  by  the  mast  increases  the  /<f ,  derivative  for  the  mastless 

vehicle  (Appendix  B)  by  270%.  It  explains  large  roll  excursions  that  occurred  in 
earlier  versions  of  the  vehicle  which  lacked  the  static  stability  of  the  current  Dorado 
vehicle. 


2 

And  note  that  the  Kr  derivative  varies  with  sm  ,  indicating  the  necessity  for  modelling 
the  effect  of  depth. 


Unsteady  (Added  Mass)  Derivatives 

The  ESAM  program  [5]  predicts  the  added  mass  coefficients  based  solely  on  input  geometry. 
Its  input  file  is  shown  in  Appendix  A.  The  specified  input  geometry  is  shown  in  black  in 
figure  6  below.  ESAM  calculates  added  mass  by  replacing  each  vehicle  component  with  an 
equivalent  ellipsoid,  shown  in  red  in  the  figure.  The  effect  of  the  velocity  field  from  the  hull 
on  appendages  is  accounted  for  by  using  hull  replacement  ellipsoids  with  geometry 
optimized  to  each  appendage;  portions  of  these  locally  optimized  ellipsoids  are  shown  in 
blue. 
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Fig.  6:  Dorado  ESAM  Input  (black)  and  Ellipsoid  (red)  Geometry 


x 


x 


ESAM  processes  the  vehicle  component  by  component.  Unlike  DERIVS,  however,  the 
mast  can  be  included  in  the  ESAM  calculation.  Although  the  mast  is  covered  with 
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self-aligning  fairings,  its  added  mass  contribution  can  be  estimated  by  assuming  the 
fairings  always  align  with  the  local  flow.  And  since  the  added  mass  of  a  cylinder  is 
primarily  determined  by  its  dimension  normal  to  the  local  flow,  it  suffices  to  add  a 
cylinder  to  the  top  of  the  hull  with  a  diameter  equal  to  the  maximum  thickness  of  the 
fairings  and  a  length  equal  to  the  average  depth  of  the  water  over  the  hull  (we  assume 
Zq  =  3  m  -  see  figure  4). 


The  propeller  guard  (essentially  a  ring  wing)  makes  a  very  small  contribution  to  the  added 
mass.  It  is  modelled  here  as  a  box  with  ellipsoids  making  up  each  side.  The  length  of  each 
side  is  1/4  the  circumference  of  the  propeller  guard. 

The  overall  added  mass  coefficients  predicted  by  ESAM  are: 


u 

-1 . 391E-03 


v 

0 . OOOE+OO 
-3 . 587E-02 


w 

2 . 187E-05 
0 . OOOE+OO 
-2 . 905E-02 


P 

0 . OOOE+OO 
9 . 979E-04 
0 . OOOE+OO 
-1 . 558E-04 


q 

1 . 957E-05 
0 . OOOE+OO 
-5.780E-04 
0 . OOOE+OO 
-1 . 822E-03 


r 

.OOOE+OO  X 
. 051E-04  Y 
•OOOE+OO  Z 
•667E-05  K 
•OOOE+OO  M 
•578E-03  N 


0 

7 

0 

-2 

0 

-1 


This  matrix  gives  the  unsteady  derivatives  by  matching  a  force  on  the  right  (one  of  the  6 
rows)  with  an  acceleration  on  top  (one  of  the  6  columns).  The  matrix  is  symmetric  so  it  is 
redundant  to  show  derivatives  below  the  diagonal.  The  component  contributions  that  make 
up  these  totals  are  shown  on  the  next  page. 


There  remains  the  issue  of  the  effect  of  variable  depth  on  added  mass.  That  is,  what  is  the 
effect  of  variable  wetted  mast  length  on  these  coefficients,  an  effect  that  is  accounted  for  in 
the  steady  state  coefficients?  This  is  ascertained  by  increasing  the  length  of  the  mast 
(effectively  the  wetted  mass  length  sj  by  10%  in  the  ESAM  input  file  and  rerunning  the 

code.  The  coefficient  changes  that  result  are  all  less  than  2%  except  for  the  Xqdot  =  Mud,n 
derivatives  which  change  by  35%.  However,  in  the  equations  of  motion,  the  Xqdot,  Mudot 
coefficients  are  always  paired  with  a  m  zG  term,  where  m  is  vehicle  mass  and  zG  is  the 
vehicle  center  of  mass  location  along  the  z  axis.  Since  Dorado’s  m  zG  term  is  more  than 
100  times  larger  than  the  change  in  X  dof  Mudot  obtained  by  lengthening  the  mast  10%,  we 
can  neglect  this  change. 

Thus,  the  equations  of  motion  are  insensitive  to  changes  in  the  added  mass  coefficients  that 
result  from  varying  wetted  mast  length.  As  a  result,  the  added  mass  coefficients  do  not 
need  to  account  for  varying  depth.  The  coefficients  do,  of  course,  account  for  the  nominal 
mast  length,  assuming  zQ  =  3  m. 
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ESAM  (version  1.42)  output  file  on  02-0ct-02  at  15:50:44  hrs . 
The  input  file  is:  ~ /MCM/RMV/Dorado/Derivs/Build3/esam . inp 
A  13  component  body. 


THE  DIMENSIONLESS  ADDED  MASSES;  the  reference  length  =  8.200000 


component 

type 

Xud 

Xvd, Yud 

Xwd, Zud 

Xpd, Kud 

Xqd, Mud 

Xrd, Nud 

Nrd 

Hull 

hp 

-0.0009087353 

0.0000000000 

0.0000000000 

0.0000000000 

-0.0000008232 

0.0000000000 

-0.0013951687 

StbdForepln 

Is 

-0.0000105089 

-0.0000001869 

0.0000037531 

0.0000002964 

-0.0000008402 

0.0000007129 

-0.0000004653 

PortForepln 

Is 

-0.0000105089 

0.0000001869 

0.0000037531 

-0.0000002964 

-0.0000008402 

-0.0000007129 

-0.0000004653 

StbdSternpln 

Is 

-0.0000162419 

-0.0000034912 

0.0000087788 

0.0000006386 

0.0000036844 

0.0000026462 

-0.0000017076 

PortSternpln 

Is 

-0.0000162419 

0.0000034912 

0.0000087788 

-0.0000006386 

0.0000036844 

-0.0000026462 

-0.0000017076 

Rudder 

Is 

-0.0000070666 

0.0000000000 

-0.0000016487 

0.0000000000 

-0.0000010920 

0.0000000000 

-0.0000338170 

VertStabilzr 

Is 

-0.0000275018 

0.0000000000 

0.0000037440 

0.0000000000 

0.0000036360 

0.0000000000 

-0 . 0001017704 

Keel 

Is 

-0.0002181763 

0.0000000000 

-0.0000038160 

0.0000000000 

-0.0000236842 

0.0000000000 

-0.0000336882 

PropGuards 

3n 

-0.0000000092 

0.0000000000 

0.0000000000 

0.0000000000 

-0.0000000002 

0.0000000004 

-0.0000042532 

PropGuardp 

3n 

-0.0000000092 

0.0000000000 

0.0000000000 

0.0000000000 

-0.0000000002 

-0.0000000004 

-0.0000042532 

PropGuardt 

3n 

-0.0000000092 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000002 

0.0000000000 

-0.0000000001 

PropGuardb 

3n 

-0.0000000092 

0.0000000000 

0.0000000000 

0.0000000000 

-0.0000000006 

0.0000000000 

-0.0000000001 

Mast 

Is 

-0.0001758907 

0.0000000000 

-0 . 0000014721 

0.0000000000 

0.0000358501 

0.0000000000 

-0.0000009597 

TOTALS 

-0.0013909092 

0.0000000000 

0 . 0000218711 

0.0000000000 

0.0000195741 

0.0000000000 

-0.0015782564 

component 

type 

Yvd 

Ywd, Zvd 

Ypd, Kvd 

Yqd, Mvd 

Yrd, Nvd 

Mqd 

Mrd, Nqd 

Hull 

hp 

-0 . 0253417740 

0.0000000000 

0 . 0000229571 

0.0000000000 

-0.0000256869 

-0.0013951694 

0.0000000000 

StbdForepln 

Is 

-0.0000129153 

0.0000886604 

0.0000071426 

-0.0000161013 

-0.0000023042 

-0.0000240652 

-0.0000028089 

PortForepln 

Is 

-0.0000129153 

-0.0000886604 

0 . 0000071426 

0.0000161013 

-0.0000023042 

-0.0000240652 

0.0000028089 

StbdSternpln 

Is 

-0.0000078196 

-0.0000516772 

-0.0000038172 

-0 . 0000219443 

0.0000035660 

-0 . 0001781714 

0.0000088000 

PortSternpln 

Is 

-0.0000078196 

0.0000516772 

-0.0000038172 

0 . 0000219443 

0.0000035660 

-0 . 0001781714 

-0.0000088000 

Rudder 

Is 

-0.0001973918 

0.0000000000 

0.0000119829 

0.0000000000 

0.0000816946 

-0.0000004150 

0.0000000000 

VertStabilzr 

Is 

-0.0005982148 

0.0000000000 

-0.0000480363 

0.0000000000 

0.0002466968 

-0.0000012560 

0.0000000000 

Keel 

Is 

-0 . 0094579572 

0.0000000000 

0 . 0010444720 

0.0000000000 

0.0003980295 

-0.0000035727 

0.0000000000 

PropGuards 

3n 

-0.0000146385 

0.0000000000 

0.0000003213 

0.0000000000 

0.0000078905 

-0.0000000001 

0.0000000000 

PropGuardp 

3n 

-0.0000146385 

0.0000000000 

0.0000003213 

0.0000000000 

0.0000078905 

-0.0000000001 

0.0000000000 

PropGuardt 

3n 

-0.0000000002 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000001 

-0.0000042532 

0.0000000000 

PropGuardb 

3n 

-0.0000000002 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000001 

-0.0000042533 

0.0000000000 

Mast 

Is 

-0.0002018634 

0.0000000000 

-0 . 0000407212 

0.0000000000 

-0.0000139188 

-0.0000085225 

0.0000000000 

TOTALS 

-0.0358679483 

0.0000000000 

0.0009979480 

0.0000000000 

0.0007051201 

-0.0018219155 

0.0000000000 

component 

type 

Zwd 

Zpd, Kwd 

Zqd, Mwd 

Zrd, Nwd 

Kpd 

Kqd, Mpd 

Krd, Npd 

Hull 

hp 

-0 . 0253417740 

0.0000000000 

0.0000256869 

0.0000000000 

-0.0000000208 

0.0000000000 

0.0000000233 

StbdForepln 

Is 

-0 . 0007270707 

-0.0000582868 

0.0001320909 

0.0000155290 

-0.0000048221 

0.0000105739 

0.0000012517 

PortForepln 

Is 

-0.0007270707 

0.0000582868 

0.0001320909 

-0.0000155290 

-0.0000048221 

-0.0000105739 

0.0000012517 

StbdSternpln 

Is 

-0.0009874937 

-0.0000733788 

-0 . 0004193077 

0.0000207162 

-0.0000058846 

-0.0000312138 

0.0000015305 

PortSternpln 

Is 

-0.0009874937 

0.0000733788 

-0 . 0004193077 

-0.0000207162 

-0.0000058846 

0.0000312138 

0.0000015305 

Rudder 

Is 

-0.0000017884 

0.0000000000 

-0.0000008424 

0.0000000000 

-0 . 0000007412 

0.0000000000 

-0.0000049594 

VertStabilzr 

Is 

-0.0000048025 

0.0000000000 

-0.0000023001 

0.0000000000 

-0.0000041904 

0.0000000000 

0.0000198096 

Keel 

Is 

-0.0002423604 

0.0000000000 

-0.0000106616 

0.0000000000 

-0.0001200035 

0.0000000000 

-0.0000439556 

PropGuards 

3n 

-0.0000000002 

0.0000000000 

-0.0000000001 

0.0000000000 

-0.0000000163 

0.0000000000 

-0.0000001732 

PropGuardp 

3n 

-0.0000000002 

0.0000000000 

-0.0000000001 

0.0000000000 

-0.0000000163 

0.0000000000 

-0.0000001732 

PropGuardt 

3n 

-0.0000146385 

0.0000000000 

-0.0000078905 

0.0000000000 

-0.0000000092 

0.0000000000 

0.0000000000 

PropGuardb 

3n 

-0.0000146385 

0.0000000000 

-0.0000078905 

0.0000000000 

-0.0000000092 

0.0000000000 

0.0000000000 

Mast 

Is 

-0.0000011091 

0.0000000000 

0.0000003770 

0.0000000000 

-0 . 0000094247 

0.0000000000 

-0.0000028078 

TOTALS 

-0.0290502406 

0.0000000000 

-0.0005779551 

0.0000000000 

-0.0001558449 

0.0000000000 

-0.0000266719 

Concluding  Remarks 

The  hydrodynamic  forces  on  the  Dorado  Build  3  vehicle  are  estimated  and  provided  in  a 
form  that  can  be  used  in  a  six  degree-of-freedom  vehicle  simulation.  The  forces  apply  to 
either  old  or  new  versions  of  the  fore  and  stemplanes.  The  drag  and  lift  of  the  snorkel  mast 
are  accounted  for  as  functions  of  vehicle  depth. 

The  force  estimates  for  all  but  the  mast  effects  are  provided  as  hydrodynamic  coefficients 
(first  and  second  order  derivatives).  The  steady  state  coefficients  are  estimated  using  the 
DERIVS/DSSP20  programs,  but  corrections  must  be  made  to  the  depth,  pitch,  yaw,  and  roll 
control  derivatives.  The  ESAM  program  is  used  to  estimate  unsteady  (added  mass)  effects. 
The  lift  and  drag  of  the  mast  self-aligning  fairings  are  estimated  using  a  rational  model.  Of 
particular  interest  is  the  lift  developed  by  the  self-aligning  fairings  in  turns  because  of 
curvature  in  the  onset  streamlines;  this  lift  is  reversed  and  amplified  by  the  self-alignment 
mechanism. 

In  calculating  the  roll  control  derivatives,  it  became  evident  that  implementing  roll  control 
via  sternplanes  is  inefficient.  An  aileron  type  control  surface  on  the  mast  would  have  more 
leverage  with  much  less  interference,  and  should  therefore  provide  better  roll  control. 
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Appendix  A:  The  Dorado  DSSP20  and  ESAM  Input  Files 


The  DSSP20  and  ESAM  input  files  listed  below  are  used  by  these  programs  for  estimating 
the  steady  and  unsteady  6  DOF  hydrodynamic  characteristics  of  a  vehicle.  This  is  done  by 
evaluating  the  vehicle  components  one  at  a  time  and,  where  possible,  accounting  for 
interference  between  components. 

DSSP20 

The  DSSP20  input  file  is  described  in  detail  by  Mackay  [4].  This  input  does  not  account  for 
the  presence  of  the  mast. 


Title  DORADO  vehicle.  Build  3  configuration. 

Plot 

Maple 

Reference 

-3.7  0.0  0.0  !  reference  is  on  hull  axis  opposite  nominal  CB 

0.0  0.0  0.0 


Hull 

Label  #Hull 
Station 


0.000 

0.000 

0.0000 

0 

1.115 

1.115 

0.9764 

- 

1 .1684 

1.1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1.1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1 .1684 

1 .1684 

1 . 0722 

- 

1.086 

1.086 

0.9263 

- 

0.94 

0.94 

0.6940 

- 

0.71 

0.71 

0.3959 

- 

0.36 

0.36 

0.1018 

0.12 

0.12 

0.0113 

0 

Nose  0.0 

Tail  -8.2  0.0  0.18 


0 

1 

I 

z  camber  relative  to  x-axis 
0.0 

009 

! 

0.0 

018 

1 

0.0 

027 

! 

0.0 

036 

I 

0.0 

045 

I 

0.0 

054 

! 

0.0 

063 

! 

0.0 

072 

I 

0.0 

081 

! 

0.0 

090 

! 

0.0 

099 

! 

0.0 

108 

I 

0.0 

117 

I 

0.0 

126 

! 

0.0 

135 

! 

0.0 

122 

I 

0.022 

086 

! 

0.067 

049 

1 

0.113 

009 

! 

0.18 

0 

! 

0.18 

Lift 

Label  #StarboardSternplane 
Sternplane  0.39  0.5842 

TailEff  0.7772  !  Dempsey 

AllMoving  0.84 
Planf orm 


-6. 

.8901 

.4842 

0. 

.  0 

-7  . 

.  1756 

1.2685 

0. 

.  0 

-7  . 

.  6021 

1.2685 

0  . 

.  0 

-7  . 

.  6021 

.2602 

0. 

.  0 

ToC  0.15 
Save 


Lift 

Label  #PortSternplane 
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Reflected 


Lift 


Label  #StarboardForeplane 
SBowplane  0.5842 
AllMoving  1.0 
Planf orm 


-1.7508  0.5842  0.188 

-1.9999  1.2685  0.188 

-2.4264  1.2685  0.188 

-2.4264  0.5842  0.188 

ToC  0.15 
Save 


Lift 

Label  #PortForeplane 
Reflected 


Lift 

Label  #Rudder 
Rudder  0.395  0.5842 


TailEff  0. 
AllMoving 
Planf orm 

.6964  ! 

0.92 

!  Dempsey 

-7.06 

0.00 

0.52 

-7.06 

0.00 

0.9 

-7.416 

0.00 

0.9 

-7.416 

0.00 

0.45 

ToC  0.25 
HasEndP  0.125  1 


Lift 

Label  #VerticalStabilizer 
Tail  0.4  0.5842 


TailEff  0. 
Planf orm 

.768 

!  Dempsey 

-6.98 

0.00 

-0.4 

-6.98 

0.00 

-1.28 

-7 . 425 

0.00 

-1.21 

-7 . 425 

ToC  0.25 

0.00 

-0.23 

Lift 

Label  #Keel 

Sail  0.5842 

Planf orm 

-3.0 

0.0 

0.5842 

-3.0 

0.0 

1 . 66 

-5.112 

0.0 

1 . 66 

-5.112 

ToC  0.16 

0.0 

0.5842 
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ESAM 

The  ESAM  input  file  is  described  in  detail  by  Watt  [5].  This  input  does  account  for  the 
mast,  assuming  a  nominal  vehicle  depth  of  3  m. 


13  Dorado,  Build  3 


P 

-3.7  put  reference  at  nominal  CB 
0.0 
8.2 
132 


hp  Hull 


/ 

18 

-8.2 

0 . 12 

0 . 12 

0 . 18 

-7.79 

0 .36 

0.36 

0 . 18 

-7.38 

0 . 71 

0 . 71 

0 . 113 

-6.97 

0 . 94 

0 . 94 

0 .067 

-6.56 

1 . 086 

1 . 086 

0 . 022 

-6 . 1214 

1 .1684 

1 .1684 

0 . 0 

-0.5842 

1 .1684 

1 .1684 

0 . 0 

-0 . 4382 

1 . 1313 

1 . 1313 

0 . 0 

-0.3895 

1 . 1016 

1 . 1016 

0 . 0 

-0.2921 

1 . 0119 

1 . 0119 

0 . 0 

-0.1947 

0 . 8708 

0 . 8708 

0 . 0 

-0 . 1461 

0 . 7729 

0 . 7729 

0 . 0 

-0 . 1168 

0 . 7009 

0 . 7009 

0 . 0 

-0 . 0584 

0.5092 

0.5092 

0 . 0 

-0 . 0292 

0.3648 

0.3648 

0 . 0 

-0 . 012 

0.2356 

0.2356 

0 . 0 

-0 . 004 

0 . 1365 

0 . 1365 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

0 . 0 

Is  StbdForepln 


-2.4264  0.5842  0.188  0 

.  0 

0 . 6756  0 . 0 

0.4265  0.6843 

0.0  0.6843 

0 . 15 

Is  PortForepln 

-2.4264  -0.5842  0.188  180 

.  0 

0 . 6756  0 . 0 

0.4265  0.6843 

0.0  0.6843 

0 . 15 

Is  StbdSternpln 

-7.6021  0.2602  0.0  0.0 

0  . 

712  0.224  0. 

.4265  1.0083 

0  . 

.0  1.0083 

0  . 

.  15 

Is  PortSternpln 
-7.6021  -0.2602  0.0  180.0 

0  . 

712  0.224  0. 

.4265  1.0083 

0  . 

.0  1.0083 

0  . 

.  15 

Is  Rudder 

-7.416  0.0  0.45  90.0  0. 

356 

0.07  0.356 

0.45  0.00. 

.  45 

0.25 

Is  VertStabilzr 

-7.425  0.0  -0.23  -90.0  0.445  0.17  0.445  1.05  0.0  0.98  0.25 

Is  Keel 


-5.112  0.0  0.5842  90.0 

3n  PropGuards 

-8.1625  .3613  0.18  -90.0 

3n  PropGuardp 

-8.1625  -.3613  0.18  90.0 

3n  PropGuardt 

-8 .1625  0.0  -0 . 1813  0 . 0 

3n  PropGuardb 

-8.1625  0.0  0.5413  180.0 

Is  Mast 

-3.2  0.0  -0.5842  -90.0 


2.112  0.0  2.112  1.0758  0.0  1.0758  0.16 


085 

0,  .3613 

.  085 

0, - . 3613 

.085  .025 

085 

0,  .3613 

.  085 

0, - . 3613 

.085  .025 

.  085 

0,  .3613 

.  085 

0, - . 3613 

.085  .025 

.  085 

0,  .3613 

.  085 

0, - . 3613 

.085  .025 

0 . 159 

0.0  0 . 

159  2. 

42  0.02. 

42  1.0 
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Appendix  B:  The  Dorado  Mastless  Steady  State 
Derivatives 


This  is  the  output  from  DERIVS  incorporating  the  corrections  to  drag  (x*),  the  small  Z  , 
Mw,  Z ,  and  M  corrections  resulting  from  the  foreplane  A  Zw  discrepancy,  and  the 
sternplane,  rudder,  and  foreplane  control  derivative  corrections. 


CORRECTED  Dorado  MASTLESS  Dimensionless  Stability  Derivatives:  Build  3 

X*_old  =  -0.0020417 


X*, Xvv 

Xw, Xww,  Xw | w 

Xpp, Xrr 
Xq,Xqq,Xq|q 

-0.0036001 

0.0000184 

-0.0000903 

-0.0001336 

0.0690062 

0 .0810214 

0 . 0048077 
0.0093780 

-0 . 0014221 

0.0001509 

X*_< 

Yv, Yv | v 

Yp,Yp|p 

Yr , Yr | r 

-0.1530390 

0.0108496 

0.0265862 

-0 .1192487 
0.0000163 
0.0039736 

Z* , Zvv 

Zw, Zww, Zw | w 

Zpp, Zrr 
Zq,Zqq,Zq|q 

-0 . 0002779 
-0 . 1480168 
0.0000000 
-0 . 0211099 

-0.0001393 
-0.0135169 
-0.0000012 
-0 . 0012814 

-0.1068810 

-0.0022532 

Zw_< 

Zq_< 

Kv , Kv | v 

Kp , Kp | p 

Kr ,  Kr  |  r 

0.0079937 

-0.0023596 

-0.0005739 

0 . 0042778 
-0.0000016 
-0.0000564 

M* , Mvv 

Mw , Mww , Mw | w 

Mpp, Mrr 

Mq, Mqq, Mq | q 

0.0004385 

0.0081326 

-0.0000116 

-0.0165276 

0 . 0020014 
-0.0020678 
-0.0002063 
-0 . 0007012 

-0.0186506 

-0 .0010120 

Mw_< 

Mq_< 

Nv, Nv | v 

Np, Np | p 

Nr, Nr | r 

-0 .0127989 

0 . 0012276 
-0 . 0070427 

0.0190693 

-0.0000021 

-0.0016514 

CORRECTED  Dimensionless  Control  Derivatives 

Xdsds , Zds 

Mds , Mdsds , Mds | ds 

-0.0022873 

-0.0165060 

-0 . 0397794 
0.0000000 

-0.0008099 

Csu 

Xdodo, Ydo 

Kdo , Mdodo , Ndo 

-0.0010989 

0 . 0017758 

0.0000000 

0.0000000 

0.0000000 

Csd 

Xdfdf , Zdf 

Mdf , Mdfdf , Mdf | df 

-0.0316288 

0.0094788 

-0 .0442960 
-0 . 0007252 

0.0006586 

Cf  = 

Xdrdr , Ydr 

Kdr , Mdrdr , Ndr 

0.0001995 

-0.0002683 

0.0056862 

0.0000106 

-0.0023979 

Cru 

CrK 

-0.1592186 

-0.0187876 


0 .0104549 
-0 .0170091 
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Appendix  C:  Mast  Drag  and  Lift  Code 


This  appendix  provides  explicit  FORTRAN  coding  for  implementing  the  mast  drag  and  lift 
hydrodynamic  forces  derived  in  the  main  body  of  this  report. 


*  Assume  zO,  xm,  rm,  1  (in  meters) 

*  theta,  phi  (in  radians) 

*  u,  v,  w  (in  m/s) 

*  p,  q,  r  (in  radians/s)  are  all  known. 

~k 

*  Statement  functions 

Fax(z)  =  (u+q*z) *sqrt ( (u+q*z) **2  + (vrxm* (l+rm**2/z**2) -p*z)  **2) 

Flat (z)  =  (vrxm* (l+rm**2/ z**2) -p*z) *sqrt ( (u+q*z) **2+ (vrxm* (l+rm**2 
#/ z**2) -p*z) *  *2 ) 

~k 

sm  =  ( zO-xm*sin (theta) ) /cos (theta) /cos (phi ) -rm 
zs  =  -rm-sm 

zb  =  -rm 

UU2  =  u**2+v**2+w**2 
ds  =  0.1871249E-2 
db  =  0.924  9517E-4 
dpmu  =  0 . 1564563E-3/UU2 
fO  =  0 . 2439302E-1 
vrxm  =  v+r*xm 

iz2by3  =  (sm**2+3*rm*sm+3*rm**2) *sm/3 

iz21by3  =  (sm**2+3*rm*sm+6*rm**2) *sm/3 

iz3by4  =  rm**2/2+rm*sm/2+sm**2/4 

lnby2  =  sm* *2 /2+rm*sm+rm* *2 *log ( ( rm+sm) /rm) 

izby2  =  (sm+2*rm) *sm/2 

rmsm  =  rm+sm 

f 0rU2  =  f 0/UU2*r 

dsfaxu  =  ds*Fax (zs) /UU2 

dbfaxu  =  db*Fax ( zb) /UU2 

dsflatu  =  ds*Flat (zs) /UU2 

dbflatu  =  db*Flat ( zb) /UU2 

lax  =  u**2*sm-2*q*u*izby2+q**2*iz2by3+sm*vrxm**2* (12*rm**3+16*rm*s 
#m*  *2  +  3*sm*  *3+2  4  *rm*  *2  *sm) / rmsm*  *3/ 6+p*lnby2  *vrxm+p*  *2  *iz2by3/ 2 
Izax  =  -U**2*izby2+2*q*u*iz2by3-2*q**2*izby2*iz3by4-vrxm**2* (2*lnb 
#y2-2*rm*sm-sm**2) /2-2*vrxm**2*izby2*iz3by4/ rmsm**2-p*iz21by3*vrxm- 
#p**2*izby2*iz3by4 

Hat  =  izby2*u*p-iz2by3*q*p+2*u*izby2*vrxm/rmsm-vrxm*lnby2*q 
Izlat  =  -U*p*iz2by3+2*q*p*izby2*iz3by4-u*vrxm*lnby2+q*iz21by3*vrxm 
Xd  =  -dsf axu-dbf axu-dpmu* lax 
Yd  =  -dsflatu-dbf latu-dpmu*Ilat 
Zd  =  0 

Kd  =  (zs*dsflatu+zb*dbflatu+dpmu*Izlat) /I 
Md  =  - (zs*dsfaxu+zb*dbfaxu+dpmu*Izax) /I 
Nd  =  xm/l*Yd 

XI  =  f0rU2*izby2* (2*v+rmsm*p+2*r*xm) /rmsm 
Y1  =  f0rU2* (-u*sm+q*izby2) 

Z1  =  0 

K1  =  f0rU2* (-izby2*u+iz2by3*q) /I 
Ml  =  -f0rU2* (Inby2*vrxm+iz2by3*p) /I 
N1  =  xm/l*Yl 
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Appendix  D:  Dorado  Build  1  Derivatives 

This  appendix  reproduces  the  corrections  to  the  derivatives  assuming  Build  1  geometry,  as 
shown  below.  The  only  differences  between  Build  1  and  Build  3  geometry  are  the  fore  and 
sternplane  shapes.  For  Build  1,  the  original  Dolphin  foreplanes  and  the  lengthened  Dolphin 
Mk  2  stemplanes  are  used. 

Fig.  7:  Dorado  Build  1  Fore  and  Sternplane  Geometry 


Note  that  both  the  Build  1  fore  and  sternplanes  use  small  end  plates  at  their  tips  (as  does 
the  rudder  in  both  the  Build  1  and  3  vehicles). 

Sternplane  Control  Derivative  Corrections 

Sternplane  Pitch  Control 

Let  CmJ  be  the  correction  that  needs  to  be  applied  to  the  DERIVS  sternplane  pitch 
control  (5v)  derivatives: 


C 


f  l  ^ 
0 


sul 


L 

\  a  J 


correct 
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1  IX(1-X) 


Then: 


l+X  5  (i  +^)2 

b  =  2.132,  d  =  .884,  *,  =  .4146 

C  =.6826 

SU1 


Sternplane  Roll  Control 

The  correction  Csd]  to  be  applied  to  DERIVS  8  derivatives  depends  only  on  X  : 


Csdl  =  3 


WBF(A»A) 
A  (1  -X2)2 


=  .4987 

Foreplane  Control  Derivative  Corrections 

The  Dorado  foreplanes  have  small  endplates  at  their  tips  with  height  h  about  5%  of 
the  full  span  length  b.  Hoerner  and  Borst  [12]  present  data  showing  that  end  plates 
increase  the  effective  aspect  ratio  a  by  a  factor: 


for: 


So: 


0^ 


ep 


1  + 


\u  7 


1.9/7 

b 

1.112 


.5  h2 
b2 


h 

7- <2.5 
b 


(\-X2)nkepb2 


r 

/  2  ^ 

1.8  +  cos(£2)/\  / 

/  k  a 

‘  4  +4 

l 2 

l  V 

cos(i2)  ) 

Since: 


b  =  2.1504  m 
a  =  2.158 
X=  .5433 
Q  =  0. 

then  the  foreplane  contribution  to  Zwl  is: 
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A  Z  ,  =  -.03933 

Wl 

Unlike  for  the  the  new,  larger  foreplanes,  this  agrees  almost  exactly  with  the 
DERIV  S/DSSP20  prediction. 

Since: 


£s  1  1  7(1  -X) 

La  =  l+\~5  (l+Xy 


then: 


=  .6271 


V 


1  lk(l-k) 


l+X  5  (i  +  X)2 


A  Z 


wl 


=  -.02467 

which  is  50%  larger  than  the  DERIVS/DSSP20  prediction  of  Zg  .  That  is,  the 

effective  DSSP20  prediction  for  Lg/  La  is  0.42  which  is  again  unrealistically  low,  as 
was  the  case  for  the  new  foreplanes. 

Since  we  have  good  Zw  agreement,  no  corrections  are  required  to  the  Build  1  values 
of  Z  ,  M  ,  Z  ,  and  M  .  The  correction  to  Zs  is  achieved  using  the  factor: 

wr  w’  q 5  q  o^.  ® 

1 

C  =-.02467- - 

% 

DSSP20 

=  1.490 
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DSSP20  Input  File 

The  Dorado  Build  1  DSSP20  input  file  differs  from  the  Build  3  version  in  Appendix  A 
only  in  the  definitions  of  the  fore  and  sternplanes.  The  Build  1  specifications  are  shown 
below. 


Lift 

Label  #StarboardSternplane 
Sternplane  0.39  0.5842 

TailEff  0.722  !  Dempsey 

AllMoving  0 . 9 
Plant orm 


-7  . 

07 

0  , 

.442 

0  . 

.  0 

-7  . 

07 

1  . 

.066 

0  . 

.  0 

-7  . 

426 

1  . 

.066 

0  . 

.  0 

-7  . 

426 

0  , 

.335 

0  . 

.  0 

ToC 

0.25 

HasEndP 

0  , 

.1  1 

Save 


Lift 

Label  #PortSternplane 
Reflected 


Lift 

Label  #StarboardForeplane 
SBowplane  0.5842  0 
AllMoving  1 . 0 
Plant orm 

-1.91  0.5842  0.188 

-1.91  1.0752  0.188 

-2.266  1.0752  0.188 

-2.266  0.5842  0.188 

ToC  0.25 
HasEndP  0.1  1 
Save 


Lift 

Label  #PortForeplane 
Reflected 
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Steady  State  Derivatives 

This  is  the  Build  1  output  from  DERIVS  incorporating  the  corrections  to  drag  (X*)  and 
the  sternplane,  rudder,  and  foreplane  control  derivatives. 


CORRECTED  Dorado  MASTLESS  Dimensionless  Stability  Derivatives:  Build  1 


X*,Xvv 

-0 . 0034084 

0 .0690062 

Xw, Xww, Xw | w | 

0 . 0000312 

0 . 0509407 

-0 . 0014697 

Xpp, Xrr 

-0 . 0000412 

0 . 0048082 

Xq, Xqq, Xq | q | 

-0 . 0001309 

0 . 0056825 

0 . 0001513 

Yv, Yv | v | 

-0 . 1530390 

-0 .1192487 

Yp,Yp|p 

0 .0108496 

0 . 0000163 

Yr , Yr | r | 

0 . 0265862 

0 . 0039736 

Z*, Zvv 

-0 . 0002779 

-0 . 0001393 

Zw,  Zww, Zw | w | 

-0 . 0820978 

-0 . 0135169 

-0 . 0844093 

Zpp, Zrr 

0 . 0000000 

-0 . 0000012 

Zq, Zqq, Zq | q | 

-0 . 0069856 

-0 . 0014308 

-0 . 0017440 

Kv,  Kv | v | 

0 . 0079937 

0 . 0042778 

Kp,  Kp  |  p  ] 

-0 . 0017005 

-0 . 0000013 

Kr , Kr | r | 

-0 . 0005739 

-0 . 0000564 

M* , Mvv 

0 . 0004399 

0 . 0020014 

Mw , Mww , Mw | w | 

0 .0174362 

-0 . 0023369 

-0 .0163539 

Mpp, Mrr 

-0 . 0000114 

-0 . 0002063 

Mq, Mqq, Mq | q | 

-0 . 0080125 

-0 . 0006765 

-0 . 0007758 

Nv, Nv | v | 

-0 .0127989 

0 . 0190693 

Np,  Np | p | 

0 . 0012276 

-0 . 0000021 

Nr, Nr | r | 

-0 . 0070411 

-0 . 0016514 

CORRECTED  Dimensionless  Control  Derivatives 


Xdsds , Zds 

Mds ,  Mdsds , Mds | ds | 

0 . 0023478 
-0 . 0076428 

-0 . 0180358 

0 . 0000000 

-0 . 0003596 

Xdodo, Ydo 

Kdo, Mdodo, Ndo 

0 . 0012530 

0 . 0007256 

0 . 0000000 

0 . 0000000 

0 . 0000000 

Xdfdf , Zdf 

Mdf , Mdf df ,  Mdf | df  | 

-0 .0192938 

0 . 0050948 

-0 . 0246659 
-0 . 0004423 

0 . 0003834 

Xdrdr , Ydr 

Kdr , Mdrdr ,  Ndr 

0 . 0001995 
-0 . 0002683 

0 . 0056862 

0 . 0000106 

-0 . 0023979 

X*_old  =  -0.0018500 


Csul  =  0.6826 


Csdl  : 

=  0 

.498 

Cfl  = 

1 . 

4901 

Cru  = 

0  . 

6281 

CrK  = 

0  . 

5602 
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ESAM  input  file 

The  Build  1  ESAM  input  file  differs  from  the  Build  3  version  of  Appendix  A  only  in  the 
definitions  of  the  fore  and  stemplanes.  The  Build  1  specifications  are  shown  below. 


Is  StbdForepln 
-2.266  0.5842  0.188 

Is  PortForepln 
-2.266  -0.5842  0.188 
Is  StbdSternpln 
-7 . 426  0 . 335  0 . 0 

Is  PortSternpln 
-7.426  -0.335  0.0  180.0 


0.0  0.356  0.0  0.356  0.491  0.0  0.491  0.25 

180.0  0.356  0.0  0.356  0.491  0.0  0.491  0.25 

0.0  0.356  0.107  0.356  0.731  0.0  0.731  0.25 
0.356  0.107  0.356  0.731  0.0  0.731  0.25 


Unsteady  Derivatives 

The  ESAM  overall  added  mass  coefficients  for  the  Build  1  configuration,  followed  by 
the  the  planform  view  of  the  ESAM  Build  1  geometry. 


U  V 

-1 . 382E-03  0 . 000E+00 

-3 . 585E-02 


w 

5 . 510E-06 
0 . 000E+00 
-2 . 673E-02 


P 

0 . 000E+00 
9 . 942E-04 
0 . 000E+00 
-1 . 405E-04 


q 

1 . 533E-05 
0 . 000E+00 
-1 . 621E-04 
0 . 000E+00 
-1 . 541E-03 


0 . 000E+00  X 
7 . 031E-04  Y 
0 . 000E+00  Z 
-3 . 033E-05  K 
0 . 000E+00  M 
-1.576E-03  N 


Fig.  8:  Build  1  ESAM  Input  (black)  and  Ellipsoid  (red)  Geometry 


x 
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